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ABSTRACT 
OO . 

We analyze the data of Kamiokande-II, 1MB, Baksan using a parameterized descrip- 
ri I tion of the antineutrino emission, that includes an initial phase of intense luminosity. 

CL(' The luminosity curve, the average energy of />e and the astrophysical parameters of 

the model, derived by fitting the observed events (energies, times and angles) are in 
$_( ■ reasonable agreement with the generic expectations of the delayed scenario for the 

C/j I explosion. 



1. Introduction 

We begin with a rapid historical excursus, with emphasis on the issues that are relevant 
^^ ' for data analysis. 

m ■ 

i,^__^ _ • Colgate &: White 1966^ propose the paradigm for the explanation of the core collapse 

(^ ' supernovae, where neutrinos are the key agents. 

OlfVI 
^ ^ ' • Nadyozhin 1978 •' concludes a detailed calculation that demonstrates an initial phase 

of intense neutrino luminosity. 



^ 



• Bethe Sz Wilson 1985 ' suggest that the energy deposition on a scale of half a second 
can re-energize the stalled shock wave. 

• 1987: Kamiokande-Il * , 1MB ', Baksaii^ and LSEMJ observe several events in corre- 
lation with SN1987A. 

• Several authors-e.g., Bahcall 1989^ -remark that non-LSD data generically meet the 
expectations. Then the main interest shifts on the relevance of oscillations. 

• Lamb & Loredo 2002 * (LL) discuss whether SN1987A data indicate specific imprints 
of the delayed scenariqfj such as the initial phase of intense luminosity. 



^With the term 'delayed scenario' we refer here and in the following to the scenario for the explosion put 
forward by Bethe and Wilson, that incorporates the initial phase of intense neutrino luminosity discussed 
by Nadyozhin (this is also called 'standard scenario' or 'neutrino assisted explosion'). 
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Figure 1; Sketch of Ve and Ue (higher and lower) luminosity curves. The excess of v^ in the first 50 ms 
('neutronization' phase) gives a small contribution to the total number of observable events. 

• Imshennik & Ryazhskaya 2004^^ suggest a 2 stage scenario with essential role of 
rotation and possible explanation of LSD events detected 4.5 hours earlier. 

Next we summarize the present status: although there is mounting evidence that the de- 
layed scenario is correct, a conclusive proof is still missing; while most people is convinced 
that SN1987A neutrinos confirmed the general picture of the explosion, there are annoying 
doubts on whether SN1987A was a standard object or not; while we continue sharpening 
our theoretical tools, we hope intensely in a Galactic supernova event to progress in our 
understanding. 

In short, one should be aware that until the theoretical picture will be definitively 
assessed, the interpretation of the results from SN1987A will continue to contain elements 
of uncertainty. But if we want to test the expectations, we are forced to answer the question: 
What are the generic expectations for the delayed scenario ? We summarize the features that 
are relevant for our analysis: 

• In conventional detectors (scintillators and water Cherenkov) the main detection 
reaction is the inverse-beta decay of electron antineutrinos on free protons (IBD): 
I'eP — *■ ne~^. 

• The main part of the emission happens in two stages as shown in Fig. [TJ 10-20% of 
the energy is radiated in an early phase, here called accretion, that should last about 
half-a-second; 80-90% of the energy is radiated later, during the phase of neutron star 
cooling, namely in a quiet thermal phase. 

• The main reactions for energy radiation during accretion are e~p — > ni^e and e~^n — > 
pi^e, the second being the inverse of IBD. The presence around the nascent neutron 
star of an abundant amount of e~^e~ plasma ensures a large flux of u^ and z^g, that 
are key ingredients to revive the stalled shock wave. [See Janka^ for a review.] 

• During the cooling phase, neutrinos of all species (i^g, '^z^, '^r, i^e, i^/x, i^r) are radiated 
with similar luminosities; this feature is occasionally called 'equipartition'. 



In this talk, based on a work in collaboration with M.L. Costantini and A. lanni^, we 
present an analysis of SN1987A data along these lines. We will compare two models for 
neutrino emission, namely the conventional one, when the neutrino emission is described 
by a 'one-component' (cooling) model, and a 'two-component' (cooling-l-accretion) model 
which resembles more closely the expectations of the delayed scenario. Preliminary results 
have been documented iiJi^lMLLMiDji j ggg ^IscM^, in particular Sect. 3.2 there. 

2. Analysis of SN1987A Observation with One-Component Model 



2.1. The Conventional Model for Ve Emission (= Exponential Cooling) 

This is the model adopted in most of the quantitative analyses of SN1987A events (the 
notable exception, to be discussed later, being the work of LL^?). The thermal emission 
of the neutron star is parameterized by a black body model with a steadily decreasing 
temperature: 



dtdE 27r I + eii.-p[E /Tdt)] (1) 

with T,{t) = T, ■ exp[-V(4rc)] 

where the suffix c means 'cooling' and where we use for convenience the natural units 
h = c = ks = ^. The 3 parameters of the model are: 

1. the neutrinosphere radius Re that describes the intensity of the emission; 

2. the initial temperature Tc that fixes the average energy of v'e] 

3. the luminosity time scale Tc that quantifies the duration of the process. 

An isotropic emission from D = 50 kpc is assumed so that the flux is simply $c = 
l/(47ri^^) X dNc/dEdt. For obvious reasons, this model for the antineutrino emission is 
called 'exponential cooling'. 

2.2. Procedure of Analysis 

Next we describe the procedure of analysis we adopted. It is significantly more complex 
than the usually followed procedures (though we verified a posteriori none of these technical 
improvements leads to qualitative changes in the conclusion) and elaborates on the results 
obtained by LlP: 

a) Since the absolute times of Kamiokande-II and Baksan are not measured precisely, and 
the time between the first neutrino and the first event in 1MB is unknown, we use only the 
relative times of the events (ii = for any detector). 

b) Since the duration of the signal is not known a priori, we consider the data in a unified 
time window of 30 s. 



CD 



CD 
> 

LU 




SCVM test for 

a£cr.+cooUng+bkgr.= 93% 
cooling+background=23% 



5 10 15 20 25 
Time [sec] 



Figure 2; Temporal distribution of the events; the end of the first second is marked by a dashed line. 



c) We account for the measured background (Kamiokande-II and Baksan), dead-times and 
angular-bias (1MB). 

d) We account for the error in the energy measurement in all detectors. More details of the 
analysis are described in Appendices A and B. 

The expected number of signal events obtained from the differential rate 



dN 



N„4^{E,, cos e)<^,St, E,)Ucos e)UEe' "^^^ 



(2) 



dtdE(.d cos 9 dcosO ' " ' dE^.^ 

that requires to know the number of protons in the detector, the interaction (IBD) cross 
section, the I'e flux, the angular bias, the detection efficiency and the Jacobian. 



2.3. Results 

A straightforward fit yields the following values of the three astrophysical parameters 



Re = 26 km, Tc = 4.5 MeV, r^ = 3.9 s 



(3) 



while the best fit of the offset times are zero, due to the fact that the signal is bounded to 
decrease with time. Our best fit is very close to the value in Bahcall book. ' 

The difference with LL is mostly due to their treatment of efficiency-not of background, 
inclusion of angles, or better cross section. Indeed, LL include efficiency in the exponential of 
e~Z^ Mj(wit ) j-| ^^(-vv/o), we include it in both terms. A formal justification of our procedure 
is given in Appendix A. 



2.4. Should We Stop Here? 

At this point we face the question: should we be content of the generic agreement with 
the expectation and stop the analysis here? One could presume that Bahcall would have 
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Figure 3; The yellow (gray) centers surrounding the NS (dark gray) represent pictorially the individual 
reactions e^p — > nve and e n — > pi?e occurring during accretion. 



answered "yes". Indeed, in his bookj- 
statement: 



the analysis of SN1987A data is introduced by the 



Unfortunately, a "minimum" model has proven adequate to describe the sparse 
amount of data. 

The study of the 'exponential cooling' model is then commented with the sentence: 

The success of this simplified "standard" model suggests that it will be difficult 
to use the neutrino events observed from SN1987A to establish more detailed 
models. 

We wish to note that "difficult" does not mean "impossible" and, remaining aware of this 
authoritative opinion, we would like to continue the discussion begun by Lamb & Loredo: 
whether it is worthwhile to go beyond the exponential cooling model. 

Actually, observations offer us some motivation to proceed. In fact, there is a hint of 
an increased luminosity in the first second. This is evident from Fig. [2l where we show the 
temporal distribution of the observed events, fronnM^ (see also^). For comparison, we 
show also two time distributions, comprising also 7 background events (the expectations is 
5.6 in Kamiokande-II and 1 in Baksan). In the lower one, all 22 signal events belong to the 
cooling component. In the upper one, only 13 signal events belong to cooling; the remaining 
9 belong to the accretion component. In the figure it is indicated the GOF of these two 
hypotheses, evaluated with a Smirnov-Cramer-Von Mises test. 

3. Analysis of SN1987A Observation with Two- Component Model 



3.1. How to Describe the Initial Emission (= Accretion) 

The initial emission is dominated by the quasi-transparent accreting region: on top 
of the previous black body emission, we simply need to model the radiation of T>e from 



e^n — > pve- Quite directly, we use: 



dNa I .. , . .„ . E: 



dtdE 'K'^ 1 + exp[i?e/-?a(t)J 

The conceptual scheme is represented pictorially in Fig. [3l where the second astrophysical 
parameter Ma (the accreting mass) besides the temperature of the positron plasma, is 
introduced. The assumed description of how Tq and Ma evolve in time is given in Appendix 
B, but it should be not a surprise that this description brings in a third parameter, Tq, the 
duration of the process of accretion. We compare in the next page our description of Ve 
emission with the one suggested by Lamb and Loredo. 

The energy spectrum of our model, for selected values of the astrophysical parameters, 
is given in Fig. U] (fronnM^). The spectrum is 'pinched' and for the selected values of the 
astrophysical parameters it has a 'pinching' factor of about ^ 4 (note that the 'pinching' is 
an output, not an input). 
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Figure 4: Continuous curve: Ve. flux for Ala = 0.15 AIq and Ta — 2.5 MeV. Dotted curve: black body distri- 
bution with the same luminosity ('l.l x 10^'' erg/s) and average energy (li MeV), namely, with parameters 
Re = 82 km and T^ = 4.1 MeV. 



3.2. Results 

Following the same procedure adopted for the one-component model, we obtain the 
values of the astrophysical parameters: 

Lamb & Loredo Model The best fit point \^^ : 

(5) 



Re = 12 km. 


Tc = 5.5 MeV, 


Tc = 4.3 s. 


\da = 5.5 Mq, 


Ta = 1.5 MeV, 


Ta = 0.7 S. 



The big value of Ma and the small value of Ta are dictated by KII early events. The 
luminosities and average energies are given in Fig. [5l The sharp transition at i ~ Ta does 
not look very appealing, when compared with the result of a typical simulation.^^'^ 



Figure 5: LL model: antineutrino luminosity and average energy in the best fit point. 
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Pagliaroli et al. Model The best fit point i^: 



Re = 16 km, Tc 

Ma = 0.2 Mp 
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(6) 



Note tliat Ma is much smaller and at the same time Tq increased. The luminosity and 
average energy curves of Fig. [6] are much more regular and in much better agreement with 
the expectations.^ 



Figure 6: Pagliaroli et al. model: antineutrino luminosity and average energy in the best fit point. 
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One may ask: what are the main reasons of the different result? This is answered 
by Tab. [T] frorrM^, where we show the effect of the various modification on the best fit 
point and we indicate the improvement of the two-component models in comparison with 
the exponential cooling model. The two main differences of our model for antineutrino 
emission with the model due to LL are 1) that we assume a time dependent temperature 
Ta, that makes the average antineutrino energy continuous; 2) that, differently from LL, we 
postpone the occurrence of the cooling phase to the end of the accretion phase. In such a 
manner we have an approximately thermal spectrum at any time, as expected,^"'^ and we 
remove the bi-modality that characterizes the first second of emission in the LL model. It 
should be noted that these defects of the LL model have been noted already in Mirizzi and 
Raffel 



Subsequent 


Rq 


Tc 


Tc 


M, 


Ta 


Ta 


Signif. 


improvement 


[km] 


[MeV] 


¥\ 


[Mo] 


[MeV] 


¥\ 


[%] 


technical (kjll) 


12 


5.5 


4.3 


5.6 


1.5 


0.7 


99.8 


Tait) 


14 


5.0 


4.8 


0.8 


1.8 


0.7 


98.9 


time shift 


14 


4.9 


4.7 


0.1 


2.4 


0.6 


98.0 


oscihations 


16 


4.6 


4.7 


0.2 


2.4 


0.6 


98.0 



Table 1: Best- fit values of astrophysical parameters. Each line of this table is an incremental step toward the 
final improved parameterization. Last column shows the significance of the 2-component models in comparison 
(likelihood-ratio test, -1-3 d.o.f.) with the exponential cooling model. 

The two-component models are better than the exponential cooling model. As a matter 
of fact, the LL-model would seem to fare better in describing the data, but it deviates more 
strongly from the expectations of the delayed scenario than our model. Indeed, the LL 
model has less difficulty to account for the large difference between the energies detected 
by Kamiokande-II and 1MB in the first second, since it can ascribe these two datasets to 
two different but contemporaneous phases of emission. This bi-modality is forbidden by 
construction in our model. 



3.3. Errors in the Pagliaroli et al. Model 

The error on the astrophysical parameters can be determined from the datsM^: 



Rr = 16+^ km 
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We also show the values of the offset times: 
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(7) 



(8) 



We see that the limited statistics manifests itself in relatively large errors. 

The marginal distributions for accretion parameters Ma and Ta of Pg and for cooling 
parameters Re and T^ of Uf, with the complete emission model are given in Ref.^^. Other 
solutions exislM^ at Ma ~ Mq. Here we discarded them, considering that in the delayed 
scenario. Ma should be a fraction of the outer core mass Mqc ~ 0.6 Mq. 



3.4- Meaning of the Individual Events 



Given a theoretical model for antineutrino emission, we can evaluate the meaning of 
each individual event. This is given in ' in table form and also in Fig. [7] in graphical form 
for the 11 events detected in the first second and using the best fit models. We derive two 
principal conclusions: 

i) First, all 11 events are mostly due to signal and not to background, with the exception of 
the event number six of Kamiokande-II that is below the energy threshold of 7.5 MeV. The 
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A Posteriori Information on the Events of the First Second: Lamb & Loredo Model 



A Posteriori Information on the Events of the First Second: Pagtiaroli et al. Model 
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Figure 7: Events in the first second: 6 in KII, 3 in 1MB, 2 in Baksan. First bar (light blue/gray), accretion 
probability; second bar (yellow/light gray), cooling probability; third bar (black), background probability. 

background probability of this event is larger in the model by Pagliaroli et al., due to the 
fact that it happened 0.686 s after the first one, namely, in the moment when the average 
energy of antineutrinos is maximal: see Fig. [6l 

ii) Second, we see from Fig. [7] that in the LL model, the 3 events of 1MB are attributed 
to cooling, the 6 events of Kamiokande-II to accretion. This is a manifestation of the 
bi-modality of the LL model noted previously: the cooling and the accretion phases are 
assumed to be contemporaneous and with very different average energies. Thus the (low 
energy) events of Kamiokande-II are easily explained by accretion, the (high energy) ones 
of 1MB by cooling. In the other model^ the spectrum is not bi-modal by construction. 
Thus the 11 events are treated on equal footing and all of them (including IMB's) have a 
high probability to be due to accretion. 

Another manifestation of the bi-modality of the LL model is the great difference between 
the time integrated spectra of accretion and of cooling; see Fig. 2 of Ref.^^'^. This implies 
that the resulting total spectrum cannot described by a thermal distribution, even including 
a 'pinching' factor. 

4. Conclusions 



Our study confirms earlier results with the 1-component model. The refined treatment 
of background, cross section, description of the scattering angle, inclusion of Baksan data, 
etc. do not lead to important changes. 

We confirm the results of Lamb & Loredo in particular the very strong evidence for 



accretion, when their 2-coniponent model is adopted. 

We discussed an improved 2-component model, where the average energy and luminosity 
curves are constrained to be continuous, cooling follows accretion, oscillations (not very 
important a posteriori) are included. 

The best fit of Tq, Ma, etc. are close to expectations; the binding energy 2.2 x 10^^ erg 
is lower than for the 1-component model; the evidence of accretion is not as strong as for 
the Lamb &: Loredo model, but still important. 
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Appendix A: Derivation of the Likelihood Function 

We present here a derivation of the Ukehhood function that follows the one given in 
Appendix A of LL, but notations and conclusions are somewhat different. For a given 
model of the detector and of the antineutrino signal, we construct the likelihood: 

/^ TT Ti , T^ I P(d = 0) if no events are detected in the i-th. bin / , -, % 

>C = Pj where Pi = < J},, J.n . ■ i x . , ■ .1 • ., , ■ (A.l) 

^_\ 1 F(a{xi)) it one event is detected m the i-th bm 

the i=th time-bin being located in 

{i - l)6t <t<i6t (A.2) 

Of course, d = means that no event was detected in the given time bin. d{x) means that 
one event was detected in the time-bin with coordinates (namely: energy, direction and 
position) generically denoted as x. The product is on the number of time bins, N, each one 
of size 5t; T = N5t. The number of bins is much larger than the number of detected events, 
N ^ Nd. The two types of probability are calculated as follows: 

P{d = 0) ~ P{d = 0,e+ = 0,b = 0)+ 
P{d = 0,e+ = 0,6 = 1)+ 
P(d = 0,e+ = 1,6 = 0) (A.3) 

Pid{x)) ~ P(d(x),e+ = 0,6 = 1)+ 
P{d{x),e+ = 1,6 = 0) 

namely by considering the various cases, when the given time bin contains or 1 events due 
to background, 6, or due to positrons produced in the detector, e"*". This separation might 
appear cumbersome at first sight but it helps for comparing with Appendix A of LL. The 
notation e^ = 1 (6 = 1) means that one positron (one background event) was in the detector 
in the given time bin, with some value of the the coordinates; similarly, d = 1 means that 
one event was detected in the given time bin, with some value of the the coordinates. 

Now we will expand the five terms in Eq. IA.3l to order 0{6t). 

We begin with the simplest case when no event is present: 

P(d = 0, e+ = 0, 6 = 0) = 

P{d = 0, e+ = 0|6 = 0) X P(6 = 0) = 

P{d = 0|e+ = 0, 6 = 0) X P(e+ = 0) x P{b = 0) = (A.4) 

1 X exp(-S'+(t)(5t) x expi-BSt) = 

1 - 5t{S+{t) + B) 

where we introduced the positron production rate S+{t), the (time independent) background 
rate B and assumed Poisson statistics for both processes. 



The next two cases have one background event but no positron produced: 

P(d = 0,e+ = 0,6 = 1) = 

' (A.5) 

P{d{x),e+ = 0,b = l) = 6tB{x) 

This is basicahy a definition of what we mean by background: The background is a detected 
event that is not due to the signal. Note that we are not attempting to describe the nature 
of the background with this procedure, we are simply taking into account its existence. 
However, having assumed that the average background is independent from the time of 
observation, there is practical way to obtain B{x); this can be measured by the experimental 
collaborations in the instants of time when the signal is absent. 

The more complicate cases are the remaining two, when a positron is produced. Adopt- 
ing the symbol X-^- for the coordinates of the positron, we get: 



P{d = 0,e+ = 1,6 = 0) = 




Jdx+P{d = 0,e+{x+),b = 0) = 




fdx+P{d = 0,e+{x+)\b = 0) x P(6 = 0) = 




/ dx+P{d = 0|e+(x+), 6 = 0) X P(e+(x+)) 


xP{b 


/ dx+{l - r]{x+)) X S'+(t, x^)6t x 1 = 




St{S+{t)-Jdx+v{x+)S+{t,x+)) = 




5tiS+it) - Sit)) 





0) = (A.6) 



where we used P(e+(a;_|_)) = S+{t,x+)5t exp(— 5+(t,x+)(5t) and introduced the definition 
of detection efficiency, that agrees with the one of LL: 

P{d = 0\e+{x+),b = 0) = l-r]{x+) (A.7) 

The last term to evaluate is: 

P{d{x),e+ = 1,6 = 0) = 

/ dx+P{d{x),e+{x+),b = 0) = 

/ dx+P{d{x), e+{x+)\b = 0) X P(6 = 0) = 

/ dx+P{d{x)\e+{x+),b = 0) X P(e+(x+)) x P(6 = 0) = (A.8) 

J dx+p{x,x+)rj{x^) X S+{t,x+)5t x 1 = 

6t J dx+p{x,x+)rj{x+)S-i-{t,x+)) = 

6tS(t,x) 

where the last equality is simply the definition of S{t,x). 
We introduced the new function p{x, y) according to: 

P((i(a;)|e+(x+),6 = 0) = p{x,x+)ri{x+) (A.9) 

The interpretation of this function is as follows. When we consider the probability that a 
positron is either detected or missed, we should get 1: 

P{d = 0|e+(x+), b = 0) + P{d = l|e+(x+), 6 = 0) = 1 (A.IO) 



In formulae and with the definitions introduced above: 



1 — r]{x^) + / dxp{x, x+)7?(x+) = 1 (A. 11) 



that imphes J dxp{x,x^) = 1. Thus, p{x,x^) can be thought of as the probabihty that 
the event produced with coordinates x+ and eventually detected, will be observed with 
coordinates x: p{x,x+) = P{d{x)\e^(x+) , e+ detected). This function is occasionally 
called 'smearing', since a positron with well-defined coordinates will be mapped by the 
detector response into a wide set of coordinates of the detected event. This concludes the 
calculation of the five terms in Eq. IA.3I I 

Plugging the results of Eqs. IA.41 IA.51 lAlBl and lA.Sl in Eq. IA.3l we find that the expression 
of the probabilities in the given time bin at the order 6t: 

P{d = 0) = (1 - 6t{S+{t) + 5)) + + 5t{S+{t) - S{t)) ~ 

^ g-5t(B+S(t)) 

(A.12) 
P{d{x)) =5t{B{x) + S{t,x))^ 

~ 5t{B{x) + 5(t,2;))e-^*(^+^W) 
These can be introduced in the formula for the likelihood finding: 

C = e-Jo '^^(B+sit)) -Q ^^(^^^^^ ^ g^^^^ ^^^^ (^_^3) 

We can omit the multiplying, constant factors exp(— i?T) and [St)^'' without affecting the 
parameter estimations. We note finally that, owing to the definitions of S{t) and S{t,x), 
it is possible and occasionally convenient to define a detector-weighted interaction rate by 
including the efficiency r]{x+) in the positron production rate: 



R{t, x+) = S+{t, x+)r]{x+) (A.14) 



The likelihood becomes 



C = const. X e-Jo dtdxR(t,x) Y[{B{xi) + / dxp{xi,x)R{ti,x)) (A.15) 



i=l 



''We remark here the main differences with LL: 
The first one emerges when we compare the equation for P{d{x)\e'^ {x+) ,b — 0) with their Eqs. (A. 20) and 
(3.21) of LL. One concludes that the function £i defined there corresponds to p{x,x+) but the factor 77(3;+) 
has been replaced with a term proportional to the step function Q{e — eo). 

The second difference regards the meaning of background. For instance, their eq. (A.14) is given by the 
product of "background" and "efficiency"; this implies that some "background" events will be seen, some 
other will be missed. This conceptual scheme is appealing being completely analogous to the one adopted 
for the signal events, but in order to be practical, it would need a model for the "background" (namely, 
contamination of neutrons, of radon, electronic noise, etc) along with a description of how the detector 
responds to these events (in particular, what are the "efficiencies" and what are the "smearing" functions). 
We believe that it is much simpler to consider as background event (without quotation marks) an event that 
we detect in absence of signals and this is the essence of the definition we described previously. 



Appendix B: Details of the Pagliaroli et al. Model 

We estimate the theoretical parameters by the x^'- 

d=k,i,b 

where Cd is the likelihood of any detector (/c, i, b are shorthands for Kamiokande-11, 1MB, 
Baksan). The likelihood of each of the 3 detectors is: 

Cd = e-Z^/^W'^* X nil^ e«(*')-'' x [B {E,, cos O^i) + J dEp{E, Ei)R{U,E, cos ^i)] . (B.2) 

Setting /rf = 1 and t^ = (that is adequate for Kamiokande-11 and Baksan) and noting that 
the IBD reaction is only weakly directional (that implies that the smearing on the angle 
can be in good approximation neglected) we see that this is equivalent to the likelihood 
discussed in detail just above. 

The time dependent temperature is described as: 

TJt) = T. + (T, - T,) (l)"'w,.h { ^;:^"^ ^^ (B.3) 

with m = 1 — 2. The drop of the number of neutrons with time is: 

Njt) = ^xMax(^yx ^-^ (B.4) 

^' mn \Ta{t)J l + t/0.5s ^ ^ 

where Yn = 0.6 and 

j{t) = exp[-(t/Ta)2] (B.5) 

The time shift described as follows: 

$P, (t) = <5„(t) + (1 - i(t)) X $,(t - Ta) (B.6) 

Finally we give some details of the treatment of oscillations: For normal mass hierarchy 
(used in the calculation given in the main text) the survival probability and the observed 
z/g flux are: 

For inverted mass hierarchy u — v interaction introduces a swap between the v^ and i>x so 

that 

P = U^Pf + U^^{1 - Pf), 

n n (B.8) 

<^,^=P<^l + il-P)<^l. 

The earth matter effect is described by the PREM model. 



Appendix C: Discussion at the Workshop 

We report the discussion at the Workshop after the talk was dehvered. This discussion 
emphasizes certain important points and it is (at least in our view) useful, but since we 
were not able to reproduce the exact words used we apologize in advance (and take the full 
responsibility) for possible misunderstanding. 

Q [P. Vogel] : Don't angular distributions disagree with expectations? In particular, 
IMB's? If you find the GOF is had, you could doubt that such an analysis is legitimate. 

A : The GOF of the angular distribution is discussed in our previous publication, PRD70:043006,2004. 
It is better than the conventional 5%, so the problem is not so severe as one may feel. Our re- 
evaluation of GOF and also the analysis uses the newly calculated cross section (PLB564:42-54,2003) 
and the published angular bias of IMB^ ' . 

Again on 1MB, there is an important point to note: the events of 1MB cannot be elastic scattering 
events, because they are not so directional. And even being ready to consider something exotic tak- 
ing place, it is hard to imagine a reaction that is forward peaked but too much, as needed to locate 
half of the events of 1MB in the region where they are, 30° < 6 < 60°. All in all, the hypothesis 
that they are inverse beta decay events seems to be the most reasonable and conservative. 
Other Authors explore the hypothesis that 1MB analysis could be biased, e.g., Malgin in Nuovo 
Cim. 021:317-329, 1998 (that incidentally also proposes interesting remarks on the time distribution 
of the events): this is not the case of our work on SN1987A events. 

Q [V. Berezinsky] : The role of rotation is important and should be included in the 
analysis. 

A : We tried to avoid linking our analysis to a specific model, and we kept the parameters free. 
These results could be used as follows: Suppose that a certain model with rotation produces an 
initial temperature of accretion Tq much lower than the one of our best fit value; then, one could 
be entitled to conclude that this model with rotation is disfavored by the observations of SNI987A 
neutrino events. 

Q [V. Matveev] : Could some events be due to neutrino interactions with iron? 

A : Yes the most complete analysis should take into account all the interactions happening in 
the detectors. We tried to keep this analysis as simple as possible and this is the reason why we 
mainly focused on the parameterization of the most important flux, namely, the one of electron 
antineutrinos. As emphasized in the paper of Imshennik and RyazhskaycM^ , in order to assess the 
importance of the interactions with iron, we would need to describe the electron neutrino flux, too. 
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